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Abstract 

Rigid particles pack into structures, such as sand dunes on the beach, whose overall 
stability is determined by the average number of contacts between particles. However, 
when packing spatially extended objects with flexible shapes, additional concepts must 
be invoked to understand the stability of the resulting structure. Here we study the 
disordered packing of chains constructed out of flexibly-connected hard spheres. Using 
X-ray tomography, we find long chains pack into a low-density structure whose mechan- 
ical rigidity is mainly provided by the backbone. On compaction, randomly-oriented, 
semi-rigid loops form along the chain, and the packing of chains can be understood 
as the jamming of these elements. Finally we uncover close similarities between the 
packing of chains and the glass transition in polymers. 

It is an enduring puzzle why a boxful of ball bearings, even when compacted by many 
taps, never packs denser than ~ 0.64 [H [2l El SI El IS] - This is much less dense than Kepler's 
stacking of triangular layers, known to every grocer as the optimal packing of oranges 
(density ~ 0.74). Yet despite their sub-optimal density and lack of periodic order, jammed 
packings are rigid and resist shear. Replacing spheres with less symmetric objects, e.g., rods 
or ellipsoids, introduces new degrees of freedom that alters the packing structure and creates 
new modes of response |HJ |9j [TU1 [TTJ [12] . Here, we describe the jammed packing of flexible 
granular chains. Although density and coordination number decrease dramatically with 
increasing chain length, the packings remain rigid. Using X-ray tomography to visualize 



1 The James Franck Institute and Department of Physics, The University of Chicago, Chicago, IL 60637 
2 The Department of Geophysical Sciences and Center for Advanced Radiation Sources, The University 

of Chicago, Chicago, IL 60637 



1 



the chain conformations, we find long floppy chains effectively partition into a collection 
of nearly-rigid elements — small loops — that then jam into a rigid packing. Randomly- 
packed spheres have often been used as a model for simple glasses [2j [13] . Building on this 
and on the fact that polymer molecules are often modeled as flexible chains |14[ [T5] , we find 
that the jammed packing of chains captures the dependence of the polymer glass transition 
on chain length, topology, and stiffness. 

Our chains are the familiar ball-chains used as window-shade pulls (Fig. [TJa, inset): the 
"monomers" are hollow metal spheres, the "bonds" are short metal rods. The monomers 
are uniform in size and uniformly distributed along the chain; they are also free to rotate 
about the backbone so that the chain does not support torsion. These ball-chains are not 
arbitrarily flexible: there is a maximum bond flex angle # max . We use the minimum loop 
size £ = 27r/^ max as a simple measure of chain stiffness; this is the minimum length of 
chain that can close to form a ring. We used an aluminum chain with monomer diameter 
a = 2.4 mm and £ = 7.5, and a brass chain with a = 1.9 mm and £ = 11.2. For each, we 
studied the packing of both linear chains with free ends, and of cyclic chains whose ends 
are connected together in a loop. 

The chains are loaded into a cylindrical cell mounted on an electro-mechanical shaker: 
long chains are rapidly unspooled into the cell end-first, short chains are poured-in few at a 
time. The chains are compacted by giving the cell discrete vertical "taps" : a single period 
of 30 Hz sinusoidal vibration with peak-to-peak acceleration of 8g, where g = 9.8 m/s 2 . 
After each tap, the pack height is measured to find the packing density p. The compaction 
process found here is similar to that of hard spheres |6| : p increases in a logarithmic fashion 
with the number of taps, and is slow to approach a steady state; we chose 10 4 taps as an 
arbitrary stopping point to measure the "final" density /?f. The packing produced is always 
disordered with no signs of chain crystallization. We repeated these measurements using 
cells of two different diameters, 4.75 cm and 2.5 cm; this had no strong effect on our results. 

Fig.[T^, plots pf versus M, the number of monomers per chain. For linear chains, pi falls 
monotonically with increasing M, from pf ~ 0.64 for M = 1 to a much-reduced asymptotic 
value /?f oo for the longest chains (M — ► do). The asymptotic density is slightly higher for the 
floppy chain (£ = 7.5 and pf j00 = 0.43) than for the stiff chain (£ = 11.2 and pf j£>0 = 0.39). 
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For linear chains, varying M changes the density of chain ends in the pack. For cyclic 
chains, end-effects are absent. As expected, long cyclic and linear chains pack at the same 
asymptotic density; however, short (but still floppy) cyclic chains (M > 2£) also packs at 
Pf oo- Still smaller (M — » £), semi-rigid cyclic chains actually pack less densely than the 
long chain limit. The trend here is opposite to that found in linear chains, where short 
chains pack denser than long ones. 

That small, floppy cyclic chains pack at /jf j00 suggests end-effects are responsible for the 
enhanced packing density of short linear chains (in the floppy regime) . We can check this 
directly by packing a mixture of linear and cyclic chains of the same length. As shown in 
Fig. []J>, pf increases linearly with the fraction of linear chains. This indicates that in a 
packing of floppy chains, ends act as non-interacting, density-enhancing defects. 

To reveal the detailed packing structure, we used X-ray tomography. Here, only the 
£ = 7.5 Al chain is sufficiently X-ray transparent to be imaged. Our X-ray source is a 37 
keV beam (GSECARS beam line) at the Argonne Advanced Photon Source. Beam size 
limits us to a relatively small sample: the chains are confined inside a cylinder of diameter 
2.5 cm (= 10.5 monomer diameters), and the imaging volume covers 5 cm (out of ~ 20 cm) 
of the packing column, containing 1000-1500 monomers. On the other hand, the high image 
resolution (27.5 //m/voxel) allows us accurately locate every monomer and trace every bond. 

Fig. [2^, shows the pair distribution function g{r) for packings of linear chains from 
M = 1 to M = 4096. The most notable change with M is in the structure of the second 
peak, associated with second-nearest neighbors. For spheres (M = 1), the second peak of 
g{r) is split into two sub-peaks, a well-known feature corresponding to two distinct particle 
arrangements [16]. The sub-peak at r = 2a corresponds to a linear trimer; for chains, this 
configuration naturally suggests three successive monomers along a backbone. The sub-peak 
at r = \ /r 3a corresponds to a rhomboid cluster of 4 particles (Fig. [2^,, inset); for chains, this 
requires contact between monomers without a shared bond. As we move to long chains, 
the r = \/?>a sub-peak quickly disappears, indicating a sharp suppression in the number of 
contacts between monomers that do not share a bond. For spheres, the mechanical stability 
of the packing is provided by pairwise contacts. For long chains, the suppression of pairwise 
contacts between monomers that do not share a bond leaves the backbone as the remaining 
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source of rigidity. But how can floppy chains be arranged into a rigid structure? 

Unlike the r = \/3o sub-peak, the sub-peak at r = 2a persists as M increases: it 
broadens and shifts slightly to r < 2a. This shift suggests a large number of bonds must be 
flexed; for M = 4096, the peak center at r = 1.89a corresponds to a flex angle of 38° (Fig. [2)3, 
inset). Measured directly, the flex angle distribution p(ff) for long chains (M = 4096) peaks 
strongly 9 max = 48°; by contrast, p{9) for short chains (M = 8) is much flatter (Fig. [2^). 
We also compute the bond-orientational correlation Cb(s — s') along a single chain 

C b (s - s') = (b(s) ■ b(s')) s ; (1) 

here s, s' are monomer labels, and b(s) is the bond that links monomers s and s + 1. Fig. [2J1 
shows C b {s — s') for compacted long chains. It decays rapidly from unity, becomes maxi- 
mum anti-correlated at s — s' ~ 6, before returning to zero and executing small amplitude 
oscillations (period ~ 10.5a). This indicates there is no long-range orientational order along 
the backbone; instead, the local chain conformations contain a prevalence of near-minimal, 
semi-rigid loops (recall £ = 7.5) (Fig. [3]). For comparison, the g(r) for M = 8 cyclic chains, 
which are rigid loops, is very similar to that for long chains save for subtle differences: 
the second peak is sharper, and is located at r = 1.85a (rather than 1.89a) correspond- 
ing to next-nearest vertices of an octagon. The g{r) of the M = 16 cyclic chain is nearly 
indistinguishable from that for long chains (Fig. [2)3). 

We suggest these strongly-flexed bonds and near-minimal loops are responsible for the 
rigidity of long-chain packings. Specifically, rigidity arises from the jamming of rigid ele- 
ments (loops) that "condensed" out of a floppy object (the chain) upon compaction. Con- 
sider a pack of minimal loops, each with £ monomers; these will be completely rigid rings. 
A ring has 6 degree of freedom; for a pack of rings to be jammed, i.e., mechanically rigid, 
there must be on average 6 independent constraints per ring. These are provided by pairwise 
contacts between rings. On jamming, each ring is on average in contact with 12 other rings; 
on a per monomer basis, each monomer will have on average 12/£ contacts, or z = 2 + 12/£ 
nearest neighbors (coordination number). Compared with frictionless spheres, where the 
mean coordination at jamming is z = 6, a jammed pack of rigid rings is much less coordi- 
nated and therefore less densely-packed. Since z ~ l/£, the packing density will be even 
lower for stiffer chains (larger £). 
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The loops in compacted chains are nearly, but not exactly, minimal in size; they are 
semi-rigid and have more degrees of freedom than rigid rings. In addition, not all bonds 
are maximally flexed, and not all monomers belong on a loop. For linear chains, the chain 
ends, being bonded to only a single partner, are less constrained than monomers in the 
middle of a backbone. Each of these effects increases the number of constraints needed 
before a pack can jam, thereby enhancing the mean coordination and packing density at 
jamming. Therefore we expect: (1) long linear chains should pack less densely than short 
linear chains, but small, (nearly) rigid cyclic chains should pack the least densely of all; (2) 
stiff chains with large loop sizes should pack less densely than floppy chains with small ^. 
Qualitatively, this is what we observed. 

In the spirit of the Jamming Phase Diagram [T7|, which links the glass transition of 
simple glass-formers with the jamming of grains, we suggest there is a similar connection 
between polymers and macroscopic chains. One well-studied aspect of the polymer glass 
transition is the variation of the glass transition temperature T g on chain length and chain 
topology. For linear polymers, as M — > oo, T g quickly asymptotes to a constant value; 
whereas T g decreases rapidly as M — ► 1. This is commonly described by the Flory form 

T g (M) = T gi00 - K/M , (2) 

where K > is polymer-specific parameter [18j. For cyclic polymers, in the long-chain 
limit, T g is the same as that of long linear polymers of the same type; but as M become 
very small, T g increases [19]. As with ball chains, the behavior of linear and cyclic polymers 
follow opposite trends as M becomes small. 

The Jamming Phase Diagram suggests a way to directly compare glassy polymers with 
jammed chains. Here, temperature T and specific packing volume v = 1/p are orthogonal 
axes; jammed/glassy states occupy the high-density, low-temperature region of the phase 
diagram, the rest are occupied by fluid states. The glass transition is the transition from 
liquid to glass along the T-axis, while jamming is the transition from unjammed to jammed 
states along the f-axis at T = jTTj. In this sense, T g (M) and Vf(M) = l//9f(M) should be 
analogous quantities. When we compare T g (M)/T gi00 for linear and cyclic polydimethyl- 
siloxane and Vf(M) /«f oo f° r t ne packing of chains, the overall similarity between them is 
striking (Fig. [2]). 
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As shown earlier, in a packing of floppy linear chains, chain ends acts as non-interacting, 
density-enhancing defects. We can then write the specific packing volume Vf in the floppy 
regime as the weighted sum of v e , contributed by ends, and v\> > v e , due to the bulk, 
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2£(v h - v e ) 
Vh M ; (3) 
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here I ~ £ is the extent that end-influences propagate into the bulk of the chain. Cast in 
this way, Vf(M) takes on a form identical to the Flory form (Eq. [2| for T g (M). Fitted to 
appropriate regimes of the pf(M) data, we find for £ = 7.5, = 2.3, £(vy> — v e ) = 2.6; taking 
£ = £ = 7.5, then v e = 1.95. For £ = 11.2, the best fit gives = 2.5, and £(v\> — v e ) = 6; 
taking £ = £ = 11.2 gives v e = 1.96. 

Finally, chains with larger loop sizes pack less densely than chains with small £. The 
analogy of Vf (M) with T g (M) suggests stiff polymers will have a higher T g than those that 
are more flexible. Indeed, for vinyl polymers with rigid side-groups that restrict chain 
flexibility, T g is higher for those with large side-groups than for those whose side-groups are 
small 120]. 

These observations do not prove the polymer glass transition must be due to the jamming 
of chains, especially at finite temperatures and with realistic interactions. Indeed, the role 
of temperature points to important distinctions between polymers and ball-chains: (1) the 
stiffness of a ball-chain is given solely by its construction, but a polymer becomes stiffer at 
lower T; (2) for the packing of macroscopic objects, thermal motion is irrelevant and the 
generalization of "temperature" is an open question [21 . Therefore the purely geometric 
jamming of chains cannot be exactly analogous to the glass transition in polymers. Some 
aspects of temperature can perhaps be simulated by varying the tapping strength used to 
compact the packing [6j[22]. Interaction can also be introduced, e.g. by making the chains 
slightly sticky with a thin coating of viscous oil. But even hard particles and no explicit 
temperature, the similarities between jammed ball-chains and glassy polymers are striking 
enough to suggest that the jamming idea captures much of the physics. Since it has long 
been a matter of debate whether glass transitions in polymers and in simple glass-formers 
are fundamentally similar, it is attractive to think jamming, whether of grains or chains, 
may provide a unifying connection. 
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We have shown long floppy chains pack into a low-density structure whose rigidity is 
chiefly provided by the backbone. This can be understood as the jamming of semi-rigid 
loops that formed when the chains were compacted. By invoking an analogy between the 
specific packing volume and the glass transition temperature, we have shown the packing of 
chains parallels the polymer glass transition in important respects. If these two phenomena 
are indeed closely connected, as our data suggest, it would be a beautiful illustration of how 
molecular geometry and symmetry, independent of the specific microscopic interactions, can 
influence the structure of condensed matter. 
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Figure 1: (a) Packing density pf of ball-chains versus chain length M after 10 4 taps, for both 
linear (open symbols) and cyclic (filled symbols) chains. Each data point is the average of 5 
trials; error bars are smaller than the size of the symbols, (b) Packing density of a mixture 
of linear and cyclic chains (both with M = 16) as a function of the fraction of linear chains. 
Solid line is a linear fit. 
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Figure 2: The pair distribution function g(r) for compacted packing of chains (£ = 7.5), 
curves are shifted for clarity: (a) Linear chains, from M = 1 to M = 4096. (b) Short 
cyclic (eye) chains, compared with linear chains of the same length, and with a long linear 
chain, (c) The distribution of bond flex angles for long (M = 4096, solid outline) and short 
(M = 8, shaded fill) chain packings; dashed vertical line indicates ^max- (d) The bond- 
orientational correlation Cb(s — s') measured along a single chain, averaged over three long 
chain packings (M = 1024,2048,4096). 
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Figure 3: Tomographic reconstruction of a packing that consists of a single long chain 
(M = 4096): (a) The full reconstructed 3D image, (b) A small slice of the packing, 
indicated by the shaded box in (a), projected onto the horizontal plane. The brighter a 
particle is, the closer it is to the center plane of the slice. Highlighted particles are arranged 
in three near-minimal loops. 
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Figure 4: The glass transition temperature T g (M) of linear and cyclic polydimethyl- 
siloxane (PDMS) (open and filled squares) compared with the specific packing volume 
Vf(M) = l/pf(M) of linear and cyclic £ = 11.2 ball-chains (open and filled circles); both 
are normalized by their asymptotic M — > oo values T gj00 and ff j00 - PDMS data taken from 
|19j . Inset shows the (T, v) plane of the Jamming Phase Diagram, with trajectories for glass 
transition and jamming. 
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